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Abstract. In the fall of 1997. during an Intensive Observation Period (IOP), the 
Atmospheric Radiation Measurement (ARM) program conducted a study of water vapor 
abundance measurement at its Southern Great Plains (SGP) site. Among a large number of 
instruments, four sun-tracking radiometers were present to measure the columnar water vapor 
(CWV). All four solar radiometers retrieve CWV by measuring total solar transmittance in the 

0. 94-pm water vapor absorption band and subtracting contributions due to Rayleigh, ozone and 
aerosol transmittances. The aerosol optical depth comparisons among the same four radiometers 
has been presented elsewhere ( Geophys . Res. Lett., 26, 17, 2725-2728, 1999). 

We have used three different methods to retrieve CWV. In a first round of comparison no 
attempt was made to standardize on the same radiative transfer model and its underlying water 
vapor spectroscopy. In the second round of comparison we used the same line-by-line code 
(which includes recently corrected H 2 0 spectroscopy) to retrieve CWV from all four sun- 
tracking radiometers. This decreased the mean CWV by 8% or 13%. The spread of 8% in the 
solar radiometer results found when using the same model is an indication of the other-than- 
model uncertainties involved in determining CWV from solar transmittance measurements with 
current instrumentation. 

OCIS codes: 010.0010, 010.1110, 010.1320, 010.7340 

1. Introduction 

Solar transmittance methods can provide water vapor abundance from direct or reflected 
sunlight measurements in spectral channels in and adjacent to water vapor absorption bands. The 
so-derived water-vapor transmittance has to be translated into columnar water vapor (CWV). 
Although this relationship is well known qualitatively, 1 it has proven difficult to quantify. 
Attempts to do so for water- vapor absorption bands in the near-infrared date back to 1912. But 
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even in the last decade there has been a steady stream ot publications on this subject. For 
example, results from ground-based retrievals of CWV using sunphotometers (SPM) have been 
reported widely (see Ingold et a/. 3 and references therein). Recently, Schmid et al 4 reported on 
CWV retrievals using an airborne sunphotometer. Instruments aboard satellites, such as SAGE II 
(Stratospheric Aerosol and Gas Experiment) and POAM II and III (Polar Ozone and Aerosol 
Measurements) use the solar occultation technique (i.e., they act like a SPM by measuring the 
solar transmittance through the limb of the atmosphere) to retrieve water vapor. 5 6 Finally, CWV 
is also retrieved from airborne (such us AVIRIS (Airborne Visible Infra Red Imaging 
Spectrometer)) and spacebome (such as POLDER (POLarization and Directionality of the 
Earth’s Reflectance) or MODIS (MODerate-resolution Imaging Spectroradiometer)) instruments 
that measure the solar radiance reflected by the Earth surface. 7,8-910 

Recent findings that the H 2 0 line intensities in the visible and near infrared portion of the 
widely used HITRAN-96 database" were in error 12 and that H 2 0 lines (especially weak ones) 
might be missing from the current databases 13,14 have sparked renewed discussion of the accurate 
conversion of measured water- vapor transmittance into amounts of water vapor. 

In the fall of 1997 the Atmospheric Radiation Measurement (ARM) program 15 conducted the 
2 nd Intensive Observation Period (IOP) to study water vapor at its Southern Great Plains (SGP) 
site. Among a large number of systems such as radiosondes, microwave radiometers, raman 
lidars. Global Positioning System receivers, and an infrared spectrometer, four sun-tracking 
radiometers were present to measure water vapor. 16 

In this paper we focus on the four sun-tracking radiometers that retrieve CWV by measuring 
solar transmittance in the 0.94-|im water vapor absorption band. The measurements were made 
between 15 September and 5 October 1997 at the SGP ARM central facility near Lamont, 
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Oklahoma (36° 36’ N, 97° 22’ W, 316 m above sea level). Dry to very humid conditions, with 
CWV ranging from 1 to 5 cm, were experienced over the three-week period. As one of the steps 
in the CWV retrievals the aerosol component must be subtracted from the total transmittance in 
the 0.94-p.m band. The aerosol optical depth (AOD) comparison among the same four 

17 

radiometers has been presented previously. 

Following the philosophy of the just-mentioned AOD-comparison we first made no attempt 
to standardize on the methods used to derive CWV from the four radiometers. We found that 
three different methods had been used in conjunction with three different radiative transfer 
models. In a second round we used the same radiative transfer model (with its underlying 
spectroscopy corrected according to Giver et al. n ) for all instruments. In this paper we will show 
the results from both comparisons. 

2. Instrumentation 

The NASA Ames Research Center deployed its six-channel Ames Airborne Tracking 
Sunphotometer (AATS-6) at the SGP central facility of ARM for this IOP. This instrument, 
described by Matsumoto et a/., 18 uses an active sun sensor to keep the instrument pointed at the 
solar disk. The central wavelengths and full widths at half maximum (FWHM) for the filters are 
given in Table 1. The Si detectors are held at a constant temperature of 45± 0.6 °C. The field-of- 
view (FOV) of AATS-6 is 4.5°. A measurement sequence was repeated every 12 seconds with 
all filters scanned nine times then averaged in the first three seconds of the 12-second period. 

At the ARM SGP central facility a CIMEL sun/sky photometer measures AOD. This 
instrument is also part of AERONET, a worldwide network of CIMEL sunphotometers [Holben 
et al., 1998]. 19 The CIMEL CE-318 points to the sun based on an ephemeris calculation and then 
fine tunes the pointing with an active sun sensor adjustment. Samples consist of triplets of 
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measurements with each member of the triplet beginning 30 seconds apart and consisting of 
eight filter measurements completed within eight seconds; the triplets are repeated at every 
quarter air mass between two to seven air masses and every 15 minutes when the air mass is less 
than two. The central wavelength and FWHM for each filter are given in Table 1. The field-of- 
view is 1.2°. The temperature of the instrument is monitored but not controlled. 

The multi-filter rotating shadowband radiometer (MFRSR) 20 has a hemispherical field-of- 
view. A band is positioned to alternately move completely out of the field-of-view and then to 
block the sun according to a solar hour angle calculation allowing a measurement of the total 
downward and diffuse downward irradiance. The difference between the two measurements is 
the direct solar component normal to the receiver, and the direct normal component is calculated 
by dividing by the cosine of the solar-zenith angle and correcting for the angular response of the 
quasi- Lambertian detector. Sampling is every 20 seconds. The central wavelength and FWHM 
for each filter are given in Table 1 . The temperature is held at 40°C. 

The rotating shadowband spectroradio meter (RSS) 21 has a Lambertian receiver and a 
shadowing sequence similar to the MFRSR; however, the detector is a 512-element photodiode 
array that receives its energy input from the focus of a prism spectrograph. Sampling is 
performed once each minute. The spectral resolution between 350 and 1050 nm diminishes from 
0.3 to 8 nm because of the prism dispersive element. The temperature is held at 40°C. 

In the following we will refer to all four instruments as sunphotometers (SPM). 

3. Methodology 

In the derivation of atmospheric transmittance, we distinguish between atmospheric window 
channels and gaseous-absorption channels. The window channels are located outside of 


5 



molecular absorption bands such as O 2 or H^O bands, and are normally used to determine the 
aerosol optical depth. 

3.1. Aerosol Optical Depth 

For atmospheric window channels the SPM output voltage, V(A), obtained when observing 
the directly transmitted solar irradiance over a small bandpass AA centered at wavelength A can 
be described by the Beer-Lambert-Bouguer attenuation law 

V(A) = V 0 (A)/T 2 exp[-mr(A)] 

where V 0 (A) is the instrument calibration constant, R is the Earth-Sun distance in astronomical 
units (AU) at the time of observation, r(A) is the spectral optical depth, and m is the relative 
optical airmass, a function of the solar zenith angle. Taking the logarithm of ( 1) leads to 

lnV(A) = ln[v r 0 (A)/?" 2 ]-mr(A) ( 2 ) 

If a series of measurements is taken over a range of airmasses m during which the optical depth 
t(A) remained constant, V 0 (A) may be determined from the ordinate intercept of a least-squares 
fit when plotting the left-hand side of (2) versus m. This procedure is commonly known as 
Langley-plot calibration. 

In (1) several attenuators contribute to t(A) : 

T(A)=T x (A)+T 3 (A)+T 2 (A)-t-T a (A) (3) 

where the subscripts R, 3, 2 and a refer to Rayleigh scattering by air molecules, absorption due to 
0 3 and 02, and attenuation due to aerosol particles, respectively. 

A refined Langley technique 22 23,24 - which uses individual airmass expressions for each 
attenuator in (3) - was used for AATS-6 but not for the other 3 instruments. The window 
channels of AATS-6 were calibrated by averaging the results of 6 successful morning Langley 
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plots performed at the Mauna Loa Observatory (MLO) in Hawaii (19° 32’ N, 155° 34’ W, 3397 

m above sea level) about two weeks before the IOP. 

Calibration of Cimel #27 (the instrument deployed at SGP during the IOP) is based on a 
transfer of the calibration from Cimel #37, the reference instrument. The intercalibrations were 
performed at Goddard Space Flight Center in Maryland on 30 August 1997 and 3 November 
1997 at midday for a period of 1-2 hours. The reference instrument itself was calibrated using the 

Langley technique at MLO in May and September 1997. 

Calibration of MFRSR and RSS was based on a robust estimate using the 20 nearest 
successful Langley plots at SGP. One of those 20 nearest successful Langley plots was obtained 
with data from the morning of 29 September 1997. A Langley plot performed with AATS-6 
during that same morning yielded calibration constants that agreed within 0.5% with the Mauna 
Loa results obtained two weeks before the IOP. This suggests that during this particular morning 
the atmosphere over SGP was sufficiently stable to yield unbiased Langley plot results to be used 
in the robust estimate of the calibration constants for MFRSR and RSS. 

Once the calibration constants V 0 (A) of the window channels are known the aerosol optical 

depth t 0 (A) can be determined from (2) - (3). The AODs obtained from each instrument were 

derived independently of one another. Although the methods to remove Rayleigh, ozone and 
nitrogen dioxide optical depths may coincide in some instances, there was no attempt at a 
uniform reduction to aerosol optical depth from total optical depth. Nevertheless, AODs (A=380- 
1020 nm) obtained during the IOP by Cimel, MFRSR and RSS agreed with AATS-6 values to 
within 0.025 (rms). The AODs in atmospheric “windows” adjacent to the 0.94-pm band agreed 

within 0.015 (rms). 17 
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3.2. Columnar Water Vapor 

The Beer-Lambert-Bouguer taw, monochromatic in its nature, may be applied over small 
bandpasses AA with negligible error as long as the spectral variation of transmittance inside the 
bandpass is small. In regions of strong spectral variation of molecular absorption, such as the 
near-infrared water- vapor absorption bands, (1) may be expressed as 

V(A) = V 0 (A )R~ 2 exp[- m (r /j (A ) + r a (A ) +r 3 (A ))]r w (X ) < ' 4 ' ) 

(Note that there is no absorption due to N0 2 in the water-vapor absorption channels used here). 
T W {X) is the band- and source-weighted water- vapor transmittance 

, . f £ 0 (A)S(A)exp[-mT B .(A)]rfA (5) 

" y j [ £ 0 (A)S(AWA 

JM u 

where r w (A) is the strongly varying water vapor absorption optical depth, £ 0 ( A) is the 
exoatmospheric solar irradiance, and S(A) is the instrument response. It should be noted that 
even if E 0 ( A) and 5(A) were effectively constant over AA, the strong spectral variation of 
Tm ,(A) is sufficient to require the band-weighted transmittance T w {l) in (4). Also, (4) does not 
follow the Beer-Lambert-Bouguer law, as 7 W (X) generally cannot be modeled by an exponential 

with a negative argument of airmass times a constant band-weighted optical depth. Hence, for 
channels in strong absorption bands, V 0 (A) can no longer be found using the traditional or 
refined Langley method. In this paper we discuss three different approaches to determine V 0 (A) 
and T W (X) in order to determine CWV from measurements in the 0.94-|xm water vapor 

absorption band. 
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Method A: Modified Lang lev plot technique 


If r^A) can be modeled by an exponential with a negative argument proportional to some 
power of the slant path absorber amount such as 


7 w (A)=exp^a(m W ) /> ] (6) 

where u is the columnar water vapor and a and b are constants, then (A) can be determined 
using a modified Langley plot technique: Substituting (6) into (4), rearranging the terms and 
taking the logarithm leads to 


In V(A) + m[r a (A) + r R (A)+r 3 (A)]= lnfy 0 (A)/? -2 ]- a(muf 


(7) 


Modified Langley plots are now constructed by plotting the left-hand side of (7) versus m b . 
Therefore, the instantaneous values of the aerosol optical depth r a (A) in the water- vapor 
absorption channels are needed. These are estimated from the SPM “window” wavelengths using 
a quadratic fit on a log-log scale of r a (A) versus A. This requires the V 0 (A) values of the 
“window” channels to be determined before constructing modified Langley-plots. It is evident 
that for the construction of modified Langley plots the columnar water-vapor amount should 


remain constant, at least for the 1.5 to 2 hour period of Langley data acquisition. 

r w (I) is typically computed according to (5) over a range of slant path water vapor amounts 
using a radiative transfer model. The constants ci and b in (6) are then found by a curve-fitting 
procedure. 3,26,27 ’ 28 Combining equations (4) and (6) the CWV is 


1 

u = — l 

m 


^ m[t R (A)+ T a (A)+ T 3 (A)] 


) J 


( 8 ) 
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In this paper we used method A to obtain CWV for AATS-6 and Cimel. For the Cimel 
instrument the standard AERONET algorithm was used: the same typical filter function (5(A) in 
(5)) was used for all instruments in the network in conjunction with LOWTRAN 7 
computations 29 to determine one set of a and b. The 940-nm channel of the reference instrument 
Cimel #37 was calibrated using the modified Langley technique at MLO in May and September 
1997. For both calibration periods the V 0 (A) values of 4 morning modified Langley plots were 
averaged. The relative standard deviations in V Q (A) were ~2%. 

For AATS-6 we used MODTRAN 3.5 vl.l 29 to determine one set of a and b values for 
MLO and several sets (covering different ranges of m w u) for SGP conditions. For 5(A) in (5) 
we used the filter function of the 941.4-nm channel as measured by the manufacturer (Ban- 
Associates Inc., Westford, MA) in February 1994. The V 0 (A) value of that channel was 
determined by averaging the results of 5 morning modified Langley plots (standard deviation 
1.2%) performed at the Mauna Loa Observatory (MLO) two weeks before the IOP. . 

Method B: Differential lamn/solar spectrum technique 

This method has been described in detail by Michalsky et al. 30 Only a brief summary is 
given here. Method B avoids the need to calibrate using the modified Langley method. Instead, it 
requires the instrument output V^(A) when viewing a calibration lamp, the lamp irradiance 
E l ( A) and the extraterrestrial solar spectrum £ 0 (A) , both convolved with the filter function 
5(A) . In order to retrieve CWV we consider the ratio of the SPM output voltages measured in 
channels in ( A [n ) and adjacent to ( A 0UI ) the 0.94-^im band 


mj = E 0 ajE L a )v L aj ^ )exp [- m ( ^ ) + Xa ^ y tR y oul y x a (x oul ))] (9) 

va o J EoiKJESKWtihJ 
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Solving for r^A) we can relate this calculated value to the radiative transfer model 
calculation of r H .(X) to derive CWV. It is important to note that in (9) most measurements and 
calculated values appear as relative values, which we can determine more accurately than 
absolute values. 2131 

We have applied method B to the IOP data obtained from MFRSR and RSS. For MFRSR 
we used the 860-nm and 938-nm channels as X out and A,„, respectively. 5(A) of all MFRSR 
channels was measured in August 1996 and again in March 1998. The 938-nm channel shifted 
towards red by 0.8 nm, but no wavelength shift was observed for the other channels. In this paper 
we have used the August 1996 measurements of 5(A). For RSS we used A OMf = 871 nm (pixel 
number 440) and A,„ = 943 nm (pixel number 458). The 5(A) of each pixel has been established 
using lasers. 21 Because method B does not depend on modified Langley plots, no parametrization 
of T w (A ) is necessary, and T w (A ) can be converted into CWV using a look-up-table. For both 

instruments we used MODTRAN 3.7 v 1.0 to create such a look-up-table of T W {X ) versus mu . 
Method C: Empirical technique 

In method C the signal in the 938-nm MFRSR-channel is calibrated for the retrieval of water 
vapor by estimating the adjusted signal, V W (X) (the signal that would be measured if water vapor 
were the only attenuator) with the MFRSR while simultaneously observing the CWV, with 
another instrument nearby. The “other instrument” is a microwave radiometer (MWR) - the 
ARM CART (Cloud and radiation test bed) MWR measuring at/=31.4 and 23.8 GHz - that 
operates continuously at SGP. 32 An empirical curve can then be formed that shows the 
relationship between V W (A) and mu. An equation fitted to this curve provides an algebraic 
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expression relating V 7 H (A) and mu so that if V W (X) — the adjusted measurement is known, 
then u can be found. 

Using the definition of adjusted signal VJX) and equation (4), we have 

V„,(A) = V'(A)exp[m(r s (A)+r„(A)+T 3 W)]= V„a)R" ! 7-..W , OW 

where V(A) is the measured voltage. VJX) is fitted by the four-parameter model shown below. 

V W {X)R 2 =V 0 (X)exp\ r a{mu) b ~ l * nu \ (ll) 

The four parameters are V Q (A), a, b , and j8, where V 0 (A) is the calibration constant for the 938- 

nm channel and a, b , and /3 describe T w (x). 

This form of T W (X) is similar to the less-complicated transmission function in (6). We tried 
the simpler form but it only worked over a small range of mu typical of dry conditions during 
the winter. To extend the applicability of the transmission function over a wider range of mu 
encompassing the entire variation in vapor over the course of a year, we were forced to add a 
path dependence term, — pmu , to the exponent b where f5 is a small correction term. With this 

addition the range of validity of (1 1) is about 28 cm > mu >0 cm. Data from 15 days of clear 
sky conditions, spanning a period from 16 January to 28 August 1997 were used to determine 
V o (A) and to develop the empirical transmission function described above. These data consisted 
of 21278 20-second samples and the parameter values were found to be n=0.541 1, 6=0.5802, and 
/?=0.003284 with u in units of precipitable centimeters. 

In method C the need to calibrate using the modified Langley method and the use of a 
radiative transfer model are both avoided. However, we have to keep in mind that because the 
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parameters in (11) are determined by comparing to the MWR, method C cannot yield an 
independent measure of CWV. 

4. Results 

In the first round, as with the AOD intercomparison work, 17 we made no attempt to 
standardize on the use of the same radiative transfer model and its underlying water vapor 
spectroscopy (required for methods A and B). In the second round we used LBLRTM 5.10 33 for 
all method- A and -B retrievals. As with the AOD intercomparison work, 17 we compared all 
CWV retrievals to the AATS-6 results. Because of the different sampling strategies and days of 
operation, this resulted in as few as 466 to as many as -19,000 samples in the comparisons. The 
results were analyzed in terms of time series (Figure 1 and Figure 2) and scatter plots (Figure 3 
and Figure 4). Statistical summaries are given in table form in Table 2 and Table 3 and are 
visualized in Figure 5. 

We found that the quality of the MFRSR retrievals (Methods B and C) deteriorates at larger 
slant-path water-vapor amounts, mu, prompting us to use MFRSR data with mu < 23 cm The 
Cimel and AATS-6 retrievals do not have that limitation and no RSS retrievals were available 
for large values of mu. 

The scatter plots (Figure 3 and Figure 4) reveal a high correlation (0.995 < r 2 < 0.999) 
among the solar transmittance methods and a somewhat smaller correlation (0.984 < r 2 < 0.986) 
with the MWR. This is because the MWR and the optical instruments, despite their collocation 
did not observe the same volume of air as the viewing direction is zenith for MWR and slant path 
to sun for the SPMs and furthermore, the field-of-view of the MWR (4.5 and 5.9 (FWHM) at 
31.4 and 23.8 GHz, respectively) is larger. Consequently, we observe that the small-scale 
variations in CWV are highly correlated among the sunphotometers whereas some of the small- 
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scale features are absent in the MWR data (Figure 2). Also, we generally find the 
correspondence between SPMs and MWR to be best around solar noon. 

In the first round of comparison 3 different models have been used for the Method-A and -B 
retrievals. The results are shown in the left columns of Figure 3 and Figure 4, and are 
summarized in Table 2 and Figure 5. The differences between each method (including MWR) 
and AATS-6 range from 2.6% to 5.3% (mas). The mean differences are within ±0.1 cm and the 
mean ratios range from 0.97 to 1.03. 

In the second round of comparison we used LBLRTM 5.10 (which includes the updated 
spectroscopy according to Giver et al n ) for all method-A and -B retrievals. Repeating the 
computation with LBLRTM 5.21 (the most recent version at the time of writing) led to identical 
results. In this second round we also deviated from the standard Cimel AERONET-CWV 
algorithm (which uses a typical 940-nm filter function for all instruments) by using the measured 
filter functions for the instruments #27 and #37. The results are shown in the right columns of 
Figure 3 and Figure 4 as scatter plots, and are summarized in Table 3 and Figure 5. 

As seen from the time series in Figure 1 and Figure 2 the changes made in the second 
comparison had a significant effect. Overall they decreased the mean CWV by 8% for AATS-6 
and by 13% for Cimel, RSS and MFRSR (method B). These decreases in CWV are consistent 
with the results reported by Ingold et a/. 3 Although we observe an even better correlation among 
the different methods (r 2 closer to unity and smaller rms differences with respect to the best-fit 
line), we now find larger biases. The differences between each method (including MWR) and 
AATS-6 now range from 3.2% to 8.3% (rms). The mean absolute differences in CWV range 
from -0.22 to 0.16 cm with mean ratios between 0.92 and 1.06. The results of method A and B 
are now 6 to 14% lower than the results of the MWR and, consequently, of method C. When 
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compared to AATS-6, the MWR and MFRSR (method C) exhibit slopes deviating considerably 
from unity. Note that the MFRSR( method C) retrievals are tied to MWR results. 

Even the spread among the results of the independent solar transmittance retrievals (methods 
A and B) has increased slightly. In terms of absolute differences (column 1 1 of Table 2 or Table 
3) the spread is now 0.22 cm (0.19 cm before) or in terms of mean ratios (column 15 of Table 2 
or Table 3) the spread increased from 6% to 8%. This shows that the result of the first 
comparison round was somewhat misleading because differences in the models obviously 
compensated for other existing biases. In other words, the remaining biases must be caused by 
errors other than model errors. For method A, those primarily include uncertainty in the 
calibration constant V 0 (A), in the filter function 5(A), in the parametrization of T w and in 
aerosol optical depth r fl (A). A detailed analysis of these uncertainties can be found in Ingold et 
al. 3 For method B, uncertainties other than model errors include uncertainty in the filter function 
5(A), lamp irradiance ratio £,(A m )/£,(A ou ,) , the relative extraterrestrial solar spectrum ratio 

E 0 (\ in )/E 0 (\ oul ), and in r a (A). A detailed discussion of these uncertainties is given by 

Michalsky et al. [2000]. In this work we have used the extraterrestrial solar spectrum compiled 
by Gueymard 34 for the method B retrievals. Using the spectrum of Kurucz 35 as contained in 
MODTRAN 3.7 increased the mean CWV by 1.3%. 

5. Conclusion 

We have in hand a large data set of CWV retrievals from four sunphotometers. We have 
used three different retrieval techniques and have also compared to a microwave radiometer 
(MWR) on which one of the techniques is based. The good agreement realized in the first round 
of comparison turns out to be fortuitous because differences in the radiative transfer models 
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obviously compensated for biases found once a single model was used for all independent 
retrievals. The spread of 0.22 cm or 8% among all independent SPM retrievals when using the 
same model is an indication of the other-than-model uncertainties involved in determining CWV 
from solar transmittance measurements with current instrumentation. These uncertainties include 
primarily uncertainties in calibration and filter or slit-function profile. The changes in HiO 
spectroscopy suggested by Giver et al. l ~ had a significant impact on the SPM retrievals, 
depending on which model was used initially they decreased the mean CWV by 8% or 13%. 
With the improved spectroscopy the CWV retrievals from the SPMs are now 6-14% lower than 
the MWR results. However, this result needs to be considered in context with all CWV 
measurements performed during the IOP. A publication showing all results from the 2 Water 
Vapor IOP is in preparation. 
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Figure 1: Time series of columnar water vapor from different instruments and methods. Top: 
initial comparison where no attempts have been made to standardize on one radiative transfer 
model. Bottom: All method-A and -B retrievals use the same radiative transfer model (LBLRTM 

5.10). 

Figure 2: Same as Figure 1 but shown are only two hours to better see small-scale variations. 

Figure 3: CWV intercomparison with respect to AATS-6. Left column: initial comparison where 
no attempts have been made to standardize on one radiative transfer model. Right column: All 
retrievals done using the same radiative transfer model (LBLRTM 5.10). 

Figure 4: CWV intercomparison with respect to AATS-6. Left column: initial MODTRAN 3.5 
was used for AATS-6. Right column: LBLRTM 5.10 was used for AATS-6. 

Figure 5: Statistics for CWV intercomparison with respect to AATS-6. Top: initial comparison 
where no attempts have been made to standardize on one radiative transfer model. Right column. 
All method-A and -B retrievals use the same radiative transfer model (LBLRTM 5. 10). 
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Table 1: Central Wavelengths (X) and Bandwidths (AX, Full Widths at Half Maximum) of 


Filtered Instruments. 


AATS-6 CIMEL MFRSR 


Xfnml AX fnm] X [nm] 

AX [nm] 

X [nm] 

AX [nm] 


340 

2 



380.1 5.0 

380 

4 






413.9 

10 

450.7 5.1 

440 

10 



525.3 5.0 

500 

10 

499.3 

10 




608.5 

10 


670 

10 

665.1 

10 

863.9 5.3 

870 

10 

859.9 

10 

941.4 5.8 

940 

10 

938.0 

10 

1020.7 5.0 

1020 

10 
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op 


Cimel 
RSS 
MFRSR 
MFRSR 
MWR CART 

AATS-6 ! 

Instrument 

g- O CD 03 > 

> 

Method 

LOWTRAN 7 
MODTRAN 3.7 
MODTRAN 3.7 
n/a 
n/a 

MODTRAN 3.5 

Model 

466 

4015 

14649 

18957 

16990 


3 

1.05 -0.03 0.997 0.06 

0.98 -0.01 0.995 0.07 

1.01 0.03 0.996 0.06 

0.98 -0.01 0.997 0.05 

1.00 -0.09 0.984 0.12 


Best Fit 

slooe inter (cm) r*2 RMS (cm) 

2.83 2.93 
3.06 2.98 
2.73 2.77 

2.86 2.78 

2.87 2.79 


Mean (cm) 
AATS-6 y 

0.10 0.08 0.13 4.4% 

-0.09 0.07 0.11 3.7% 

0.04 0.06 0.07 2.6% 

-0.08 0.06 0.10 3.4% 

-0.09 0.12 0.15 5.3% 


Differences (cm) 

mean stdev RMS RMS(%) 

1 .03 0.02 

0.97 0.02 

1.02 0.02 

0.97 0.02 

0.97 0.05 


Ratios y/AATS-6 
mean stdev 


Table 2: Comparison of CWV from different instruments and methods to AATS-6. This is the initial comparison where different 
radiative transfer models as indicated have been used. 




Table 3: Same as Table 2 but the same radiative transfer model has been used for all method- A and -B retrievals. 







Local Time [h] 

Figure 1 


Submitted to Applied Optics on June 9, 2000 










18 September 1 997 


Local Time [h] 


AATS-6 (A) 
Cimel (A) 
RSS (B) 
MFRSR (B) 
MFRSR (C) 
CART MWR 






18 September 




1997 



Local Time [h] 















CWV (cm) MWR CART CWV (cm) MFRSR (Method C) 






Figure 4 


30 








CWV difference or y-intercept [cm] CWV difference or y-intercept [cm] 



CWV ratio or CWV ratio or slopo 







